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1, INTRODUCTION 

Teleoperation systems extend the human operators' capability to perform task at the remote site. In 
general, the teleoperation systems can be categorized into bilateral teleoperation[1] and multilateral 
teleoperation[2], [3]. To properly control the multilateral teleoperation system pose some challenges, 
especially when we consider a single master with multi-slaves configuration. Generally, there is inadequate 
degree-of-control on the master side to effectively manipulate the slave side, hence requires the slaves to 
have a certain degree of self-control. To address the issues related to the multilateral teleoperation, we focus 
on two main questions: |)How can we treat a multilateral system as a bilateral system?, and 2) Is there any 
mechanism that can be implemented to improve the transparency (i1.e., to provide a faithful transmission of 
force and velocity) between the master and multiple-slaves in a 4-channels structure? 

The facts that many works have been done in the past focusing on the bilateral teleoperations[4|—[6], 
then converting the multilateral teleoperation system to suit the existing bilateral structure, offers some 
advantages in terms of the controller design. Furthermore, as the total dynamics of slave manipulators can be 
changed (i.e., adding or subtracting slave manipulators, or to establish new configuration on the 
communication topology between agents), it is certainly cost inefficient if the dynamics of the master 
manipulator has also to be redesigned whenever the dynamics of slave manipulators have changed. 
Therefore, by introducing a single tunable equalizer to the master side renders flexibility to adjust the overall 
dynamics of the master manipulator to achieve symmetry between both sides of teleoperation. By this 
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approach, the system designer requires only minimal knowledge on the dynamics of the newly-changed 
slave manipulators. 

In this paper, we focus on the passive decomposition technique[7] to convert the multilateral 
teleoperation system into the bilateral form. Since the whole teleoperation process can be separated into two 
operations, namely grasping and handling, we use the sub-dynamics from the result of decomposition to 
determine the time (or in our case, ), that differentiates the time inverval between the two processes. In 
addition, we investigate the potential of introducing an equalizer connected to the master manipulator in 
order to achieve impedance matching between both sides of teleoperation. 

Due to the orthonormal properties of the Laguerre function, the equalizer can be realized by a 
Laguerre network of high order of moderately damped system[8]. In most applications, the optimum 
structure of this network can be obtained by simply selecting the optimal values for the basis of the Laguerre 
filters[9]. The equalizer should be properly tuned to improve the overall performance. Hence, we adopted the 
so-called Fictitious-Reference-Iterative Tuning (FRIT) algorithm in our problem formulation. Unlike other 
tuning algorithms which requires mathematical model of the plant prior tuning, FRIT utilizes the measured 
input-and-output data to simultaneously attain the plant model and the optimum controller parameters[10]. 
Interestingly, it requires only a single-shot of measurement data for this purpose[11]. 


2.2 PROBLEM FORMULATION 
2.1. Notations and graph theory 

Let 8 and 8” denotes the set of real numbers and vectors of dimension 
1=[1,...,1) 


It” respectively. Let 


I, © R"*” to denote 


a 
be a vector with all elements of ! and | is the transposition. We use 

! Qa) Op) ge lay ‘ ap 

the identity matrix. To denote the inner and cross products for ©: <a RP we use UD at and U XW’ 


>— 4/9)! 9) ay{ ke 2. . 
respectively. The vector norm is defined by lvl] == Vole , and oR) Nie implies 








K 
v(k)|I3 = S> lle)? (1) 
k=1 


For A = (a;;) € R"*"™ and B € R'**, the Kronecker product is defined by A @ B = (a;;)B € R™*™*. 
For number of agents, a graph G is built upon a finite set called a vertex set mapped by 
V = {v1,v2,...,Un} and the edge set of €. We say w;,,; € € if there is a communication link between 
them. Additionally, N; = { U5 |i, Vv; € E , V7} is the neighborhood set of Ui. The relationship between 
v;,v; € V for an undirected graph G is encoded by the adjacency matrix A(G) = (a;;) € RY *’ where 


QO otherwise. 


The degree matrix D(G)<¢RY*% is a diagonal matrix that D(G) =diag(A(G)1) with 


= |Ni= yar ajj;, and |\V/;| is the cardinality of A/;. The Laplacian matrix associated with G is 
defined by [12]: 


J 
L(G) := D(G) — A(G) > 0. (3) 


2.2. Teleoperation architecture 
We consider a simple agent as the master manipulator governed by the equation of motion 


MmPm = tm fe ioe (4) 


Conversely, the slave manipulators depicted by a network of agents such that each agent 2 moves 
dynamically through 


MiDi = fi = ie WiEG. (5) 
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The variables m; € IR, and p; € R? are the agent's mass and position. The subscripts / = m,7 stand for the 
master and slave manipulators, respectively. The signals f,,, € IR? and f;, € IR* are the exerted force of the 
master manipulator and the corresponding applied force of the human operator. Similarly, the signal f; € R? 
is the exerted force of each agent 7. The force f., € IR? is the estimated environment’s force felt by each 
agent 2 at its position relative to the center-of-mass (CoM) of the environment. 

The object in the environment is considered to be a solid disk of radius that moves on a sliding 
flat surface. Assuming that all agents are friction-less (1.e., they move on friction-less wheels), the friction 
force only exists between the contact of the disk and the sliding surface. The desired final position of each 
agent 2 on the environment is defined as 


: T 
pi:=R| cosa;(t) sina,(t) |, (6) 
where a;(t) = a? + y(t) is the desired angular position of i, and a? = 27(i —1)/N € R, Vi. Meanwhile, 
(o(t) 1s the environment's nominal angular displacement used in (7). For brevity, we drop the time notation (t) 


Operator Local Hand Controller : Communication Feemeote Telerobot | 
( Aurriaar (Master ) t hia Ls : (SLawe) Remote 





(a) Impedance Matching between 
Master and Locked System. 


~ 








(b) Velocities Matching through 
1. 


Figure | General Structure of Teleoperation Indicating the Equalizer § Figure 2 Concept of Impedance 
fF and Intervenient Impedances, Ci. & Cse. Matching by FRIT. 


in most equations, but use so whenever to emphasize its time-varying property. Defining F’y as the friction 
force, the dynamical equation of the object is denoted by 


N . 1 — N - 
Mepe = >) fi -— Fy, 5 meR¢ = Sj. (7) 
i=1 : i=1 
where m, € R and p,. € R? are the mass and CoM's position of the target/object, respectively. The force 


fi © R* and torque 7; € R in (7) are defined by 


(fi, —Di) 2 
_— law. aN (—pi) if |i _ Del = it 
fi= ¥ (-Bi, -pi) 
0 otherwise 


and 


oOo © 


Di fi — FP; ell nile 
JB) [o 07] italic 


0 otherwise. 


=F 
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3. MULTILATERAL TO BILATERAL TELEOPERATION 
3.1. Total Dynamics of the Slave Agents 
Following Kawakami & Namerikawa[13] and to use (6), we introduce a change of variables as 


Di = Di — Di- (8) 


The force jf; in (5) can be separated into two control components depending on the time interval, 1.e., 
grasping and handling. We express 


ie t<Thete : (grasp) 
fi = ve f (9) 
fi +ff°, t2=Ttete : (grasp + handle) 


where fi “and f/° are the formation and teleoperation control laws, respectively. The consensus control law 


fic =; ic {—h(B — pe) — b( (p; — De) )} + . Qi; {-h (pi — Pj) — b(p; — py) (10) 


=! 


is used to achieve simultaneous grasping. Here a;, = 1 if 2 can measure its relative positions and velocity 
with respect to object's CoM, and «a;. =O otherwise. The gain @;; corresponds to the communication 
topology defined in (2), while the variables k. b € IR are the control gains. The gain +; € IR is defined by 


1 


ree) Sk "1 ij 


Substituting (8) and (10) into (5) and since f,, = Ofor t < 7}.1-, we now rewrite the dynamics of agent 2 as 
MiDi = 4; ic {—k(B: — pe) — b( (p; — De )} + » Ai; 1-H ‘(pi — p;) — b(p; _ p;)} | (12) 


Suppose the target object acts as the virtual (N+1)th agents, then we _ have 
Pe = PN+1 = PN+1 — Pn+i1 With py4y1 = 0 and aje = aj(n +1). In addition, the mass of the virtual agent is 
to take the average mass of all agents given by my44 = + , m;. By this, (12) can be written as 


N+] N+1 


— pi = + ai; (pi — Dj )-bS” ai; (bi — p;). (13) 
j=1 





Figure 3. Structure of the Laguerre Network Figure 4. Communication Topology between Agents 
(No.6 is the Virtual Agent Representing the Target 
Object) 


The total dynamic consisting of Vv + 1 agents is given by 


Mp + 0(£ @ I2)p + k(£ @ I2)p = ‘ 0,  t<Trete 


f—f., t>Trete (14) 
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A A A 2(N+1 cl te OUT . . 
where p = [p],...,By4,]' € R” +1) f= [ft rey hae | fe = [Fl ove Foe]? and 


M := diag (71/41, --+ , 40/441) @ Ip. (15) 


Meanwhile, £ is the Laplacian matrix composed of the assorted actual and virtual slave agents 
defined by £L = (l;;) € RATY* AT) where 


N+1 

“ Qik i= 4 : 

lig = ~ 16) 
aij if 2 = a, 


3.2. Passive Decomposition[7] in Handling 
Introduce a new coordinate system in the form of 


+ 


z= [z, , 29 stage” 


= Sp (17) 


where S is the transformation matrix defined by 





M, My ... Mn Mn 
I Oi: - xs 0 —]J 
Ce 0 I om 0 _J = Ree tery (18) 
0 | ere I —J 
mM; Saag? 
M,; = —+Iy, j=1,...,.N+1, mz, = a 
ML; | > Vi 


The sub-coordinate z, € R* corresponds to the Locked system's coordinate. Accordingly, 
Ze = |Z ,---,2%n4,]' € R?% represents the Shape system's coordinate[7]. Substituting (17) into (14), we 
obtain the total dynamics in the new coordinate system as 


S-™MS—1%+ 6S—'(£L@I2)S12+kS~'(L@IL2)S-'z = S—' (f-f.). (19) 
With this transformation, it is observed that the mass matrix (19) is block-diagonalized such that 


S-™s7! =: | i | 


0 M (20) 


AT — vr! 2Nx2N 
where M=M >~0€R™**” Similarly, the information graph (16) is decomposed to attain 


= 
S-'(L@Ik)s' = O2x2 D |. 


21 
On x2 L ( ) 


where the structures of D € R?’*? and L € R? *2% are identical to the one defined by Lee & Spong[7]. 
The Locked and Shape systems of the total dynamics (14) as ¢ > 7}..;, are obtained as 


mii1+bD'z.+kD' 2 = fie — fer” (22) 
Mz, + bbz. + kLz. = fsn — fSp” (23) 


“ 
y ee vac tied en vt mT | =. ore 

where if Bete Sh =5'f, andl’ *° ’ Ish ? “, Notice that the Locked system (22) is coupled 

with the Shape system (23) through D In one approach, Nam & Namerikawa[14] introduced a decoupling 

control input into the system to eliminate the coupled term. However, Lee & Spong[7] emphasized that D 


Data-Driven Impedance Matching in Multilateral Teleoperation Systems (Mohd Syakirin Ramli) 


718 Oj ISSN: 2502-4752 


can be eliminated through proper selection of graph topology (3). Following result in [7], the Locked system 
(22) became the newly decoupled Locked system defined by 




















re te NW ; 
miZ = fic—fre - (24) 
 & i UEAEOURRERSTEE! 
6 
4 
2} (2) 
2 , ee 
m0 4) a {obiect 
a | "hie 
ap} 
-4| 
6 © ° 
aS See é | J 
; — - time|s] time|s| 
(a) (b) 


Figure 5. Asymptotic Convergence of the Agents onto Their Intended Final Positions on the Target 
Object as ¢ < Ti.7-, with »(0) = 0. (a) Agent’s Trajectories, (b) Velocity and Position of the Shape 
System. 


The following lemma provides the sufficient conditions for selecting gains k, to ensure that all 
slave agents simultaneously arrive at the intended final targets by the consensus control law (10). 

Lemma 1. Suppose for the total dynamics of the slave agents in (14) and ¢ < 7},;,, there exist gains 
k,b > 0 to guarantee simultaneous asymptotic convergence of all slave agents to their intended final 
positions on the environment/target iff the following conditions are met 


2b°Lsym —kM > 0 and Lan ~ 0, (25) 


_1 
sym — 9 


2 
([+L'). | L 
where is the symmetric component of “”. 
3.3. Tele-operation time Ltele 
The estimation of ‘tele is crucial since it is considered as the waiting time by the master manipulator 
before the handling session begins. This parameter is non-unique since it depends on the positions and 
velocities of all slave agents at cM), 
Theorem 1. Suppose for any * and ® satisfying Lem.1, with any arbitrary initial positions 
2N - 2N rp 
(0)€R and velocities p(0) €R of the agents, then the approximate “maximum” time ftele taken by 
the agents to reach to the intended final positions on the environment such that for a small scalar € ~ 0. 
||| =—E for i> Tele. 
ryt po ex: env —__ . . 
Proof. As for ® < Ltele, we have fsn — Ssh = Then the Shape system (23) can be written in the 
ene |r T TVS 2N 
c=[|z.,zZ.| €R and A is 


a 


form of autonomous LTI system of a(t) = Ax(t) where a stable matrix. 


1 At Lp—at 
Hence, there exist @ > Y and ke > 0 such that|l€ | < ke . Then, 


|2x(t)|| =Ile4*x(0)|| < lleA*|I|«(0)|] < ke-™*|}x(0)]). (26) 
LL... . l.; 7 
Hence, when t > = In k||a(0)||/€, ||a(t)|| < €. Therefore, Trete = = In k||a(O)|| /€.m 
4. IMPEDANCE MATCHING BY FRIT 


Consider the time interval ¢ > 7;,.;., and based on the block diagram in Figure 1, the following 
control laws for the sub-dynamics of (24) and (23) are used: 
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Tie a Cipm ans Crit i C3 fh, (27) 
fsrn= fsn- (28) 


Rewriting the dynamics of the master and slave (Locked) manipulators in terms of their impedance yields us 
to 


Pm pA and 4 = Ay Fie (29) 


where 2m = Mins and2@Le = Mxs, respectively. An equalizer is placed in the feedback loop of the master 
manipulator. By optimally tuning this equalizer will symmetrize the impedance in both sides of teleoperation. 


Contrary to Tsuji et.al.[15], we consider a special structure of f’ and the tuning procedure is carried-out 
using the 2_dimensional space data. Figure 1 also illustrates the placement of the equalizer fF in the general 
structure of our teleoperation system. The control gains C1,C2,C3 and © can be designed accordingly to 


99 CG 


form either “position-position’, “position-force” or “optimized” teleoperation architectures[16]. To avoid 


acceleration measurement, additional controllers Cme and se (low-gain PD Control) are employed as the 
intervenient impedance to the original systems[16]. Hence, (29) is rewritten as 


PatZe in On Sv 27 fe (30) 


where Bak = Zm + Cine and Zh = Zio + Cse. 


4.1. Fictitious Reference Iterative Tuning (FRIT) 

FRIT is a simple yet powerful tuning algorithm in the sense that we require only a single-shot of 
experimental data of the closed-loop system to obtain the optimum controller parameters[10]. Figure 2 
depicts the main concepts of impedance matching using FRIT in Tsuji et.al. [15]. As illustrated in Figure 


3(a), the ultimate aim of introducing the equalizer fF is to achieve 2F —* 2Lc, However, before /” can be 
tuned, it is necessary to match the velocities of the master and the Locked system through HL as depicted in 
Figure 3(b). Hence, we solve the optimization problem: 


min di, (31) 


where for the recorded initial data py, (k) and 2/(k), 
Jur =||42(k) — H(z)Bm () lic (32) 
Based on the optimal H(z), we solve the second optimization problem: 


min Jp, (33) 


where for the recorded initial data f° (k) and f?..(k) when the initial equalizer F(z) is used, 


Jp = Il f2c(k) — frel®)|l¥; (34) 


and the fictitious reference signal is 


fick) == H(z) (u°(k) + F(z) p(k) 5 (35) 
u(k) = fin(k) — Fo(z)Pm (A). (36) 


Theorem 2. Suppose that H* :=argminy Jy and F* :=argming Jp attain Jy =O and 
Jp = O,sespectively. Then, if the signal 2!(k) is rich, 


Zp=Zy,+F*= Zr. (37) 
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4.2. Laguerre network 
We can approximate a discrete time SISO system to use a Laguerre function 





1 —az)*"! 
Oe a (38) 
(z = a)" 
as UZ) = gl2julz) = pe c;L;u(z) (Figure 3). Here a := e PT and 6 := \/(1 — a?)T. The parameter 
is determined by the system's designer and ‘J’ is the sample time. Meanwhile c; € IR, 2 = 1,...,/M are the 
coefficients that form the basis of the Laguerre network. Let /;(k) © R, i=1,...,1 Mf represents the output 


of the 2th-order filter in the Laguerre network. Then, the SISO state-space model is given by 


nee { Ay 1) = i + bu(k) (39) 


where | = [l,...,ly]° €R¥,ae€ RY*™” bER™ cER™, 


| —1)'-37 ITN (n+ 7 | 
{ae os Me OM ifi>j (40) 
0 otherwise, 


| 
| 
ells 
ese. 
— 
a 
im 
—" 
— 


I= 


cf = Gi cae ca | ; (42) 
; Gesaas | Degas | / 2; 

1 =e) mM =T+-(e 7-1), 73=-Te? —-(e PF -1), m= a er —7}). (43) 
p p p 


The MIMO transfer matrix G(z) = I, © g(z) has a state-space representation given as 


G(z) : { l(k+1) = Al(k) + Bu(k) (44) 


g(k) = C11(k), 


aa heee l= 1, ® le RM A = I, Qaec R?™! xnM n= L, @ bE RM! XN and ‘eus aa I, ® cl E RVxnM 


The appropriate design of the equalizer ! can be obtained by assigning (44) into I’ To ensure faster rate 
convergence of the Laguerre network, Wahlberg[17] suggested that the value of P in (43) should be chosen 
close to the inverse of the system dominant time constant. However in our case, the parameter of ? must be 
chosen so that the designed equalizer will be in-sync with the Shape system (23). Hence, the parameter P is 
to use the parameter’s value of “ in (26) such that ? — a Therefore, we are left with the tunable parameters 
of © where the optimum equalizer ! can be obtained by simply tuning the basis (42). 


The selection of the equalizer # to solve the optimization problem of (31), requires # to be 
inverse-able and H~' is stable. Therefore, we select a biproper transfer function in the form of 


L+a@jz°>'+---+4pz-" 


h(z) = 7 == 
1+ bz +++++bDpz 


to yield us a MIMO transfer matrix of H = I, © h(z). The estimation of the coefficients 4, b; € R in (45) 
can be obtained by minimizing (32) through least-square method (LSM). Re-arranging data such that 
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po(P) ++ p81) } -a0(P) oe 29) 
Wy _ ' a RMA-P)xnP 
Bma(K — 1)-+- Bh, (K — P)-29(K — 1)--- -23(K — P) 
27(P +1) — p? (P +1) 
Ts _ E RUA-P) 


27(K) — Bin (K) 


a “ T ; : : ; 
and XH = GQ, «+: Gp by «+++ bp € IR”, to give us a linear relation of Vyy47 = Vy. Then, 


the optimal coefficients are obtained by 


XH = arg min Ju =Vy'TH (46) 


—1 
where V Hy! = (Wu H) Wy is the Moore-Penrose pseudoinverse of W;;. Next, we illustrate the tuning 
procedure to obtain the optimum equalizer /’. The cost function (34) is simplified to 


Jp = ||Fbc— Ho (u(R) + FC) Mae = |] (EE AB(E) — w(K) - FCI] 47) 


, 


where #4* is #4 with the optimized parameters obtained by (46). The overall algorithm to obtain optimum /*’ 
is summarized as following: 


Step I: 


Step 2: 
Step 3: 


Step 4: 


Step | 


Step 2 : 
Step 3: 


Set a Laguerre network F°(z) using (63) and Co = I, ®@ co € R"*"™ with an appropriate 
initial C = Co as 


1. f Wk+1) = Al(k) + Bom (h) 
oie { §i(k) = CUM(k). 


0 


; ase -a rc 0 
Perform an experiment to generate the initial datap,,,, 29, f?,. and 


Th 
Choose a bi-proper transfer function #4, and minimize (32) by (46) using data obtained in Step 
2 
Using the optimally tuned 7{(z), find the minimizer #’*(z) to minimize the cost function (47). 
The sub-algorithm in Step 4 to tune /’ is as follows: 

Generate the new data y(k) = H*' f?.(k) — u°(k). 

Using De, in Step 2, generate the state data 1(k). 

Find the minimizer c* of the constraint minimization 


° = 2 ‘ 
min |C71(k) -- y(k)| Py ie= co ||? 
subject to C? =I, @c’ EC R™*™ (48) 


This is equivalent to 

min 3 I|c* U(k) _ ji (k)|| , + n|le— coll’. (49) 
i=1 

Then, the minimizer is obtained by 

c* =TpV¥7(UpV;) (50) 

where 


Wr =([UP+1),-- UK), UP +1), UK), UP +1), ++, UK), Sal] € RMP) 
Tr =[fi(P+1),--, Hi(K), (P+1),-+, G2(K), in(P +1), +++, Gn(K), oo | ER OPM), 
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5. NUMERICAL EXAMPLE 

A numerical example to illustrate the effectiveness of our approach is presented in this section. Let 
the slave manipulators comprised of V = 5 agents with their communication topology as shown in Figure 4. 
The dashed-arrow indicates that only agent 4 can measure the position and velocity of the target object 
(indicated by virtual agent 6). The mass of the master, slave agents and the environment were chosen as 
Mm, Mi, Me] = [2.0, 0.5, 10] kg, V2. Meanwhile, the radius of the target object/environment was chosen as 




























x-direction 
SaSee7cse fief 
© ‘ 
S) \ ; 
=e = = Sez 
an 5 10 15 20 25 30 
_ y-direction 
Z. 2005 ) —— 
2 wy ee Pm iP 
5 0! aor 
“a AMM ane i ile 7 i a i a i a 
S- 5 10 15 20 25 30 
— —_ Error Norm 
= | 50 
=e 0; 2m 5 10 15 20 25 30 
<= 5 10 5, 20 25 30 a ae 
az time [s] = time [s] 
Figure 6 Velocity Matching Through Gain Figure 7 Comparison of Optimized Fictitious Signal 
(35) with the Initial /z< 


R=1m To satisfy Lemma 1, the spring and damping constants in (12) were selected as 
[k:, 6] = [1200 N/m, 500 Ns/m], With this selection, the dominant eigenvalue of the Shape system (23) has 
been obtained as Atele = Amin = 2-4029 | For ease of analysis, we assume that the friction force,” f ~ 0. 

The initial positions of the master manipulator and the C.O.M of the target(environment) were set at 
(0.0) and being stationary at ¢ © 10, Tiete), Meanwhile, the slave agents’ positions pi(0) € R°,Vi were 
initialized around the target in the region [-10, 10} x [—10, 10] € R* with zero velocity at ' = 0 ie, 
pi(0) = 0, Vi, Tp Figure 5(a), the agents’ initial positions and trajectories towards their intended final 
positions on the object for & © 10, Trete) are illustrated. Subsequently, Figure 5(b) depicts the trajectory of 
Shape system's velocity and position. It is clearly observed that the trajectories of the Shape system were 


asymptotically converging such that 17.21. Ze = lim:+7,... Z =9 The estimated teleoperation time 
based on the given initial positions were obtained as Ttele = 4.1331 ¢° Th this example, we used the 
Y — mr { 20 Y — J. gy (9 6 ; = [5 & i 

Cme = 12 ® (2 + =) Csc=h® (20 + al the controller Cm =he (100 + 3), 


p Zn = Iz ® (0.58 + 70 + =m 7 =0.1 and Tsim = 30{sec], respectively. 


intervenient impedances 


with human impedance o 
The velocity matching performance between Hpm(k) and (kt) as presented in Figure 6 verifies 
that there exists 74 such that the cost function (32) is minimized. In Figure 7, we exhibit the comparison 


‘O (J. 
between / Lc(k) and the response of fictitious signal (35) when the equalizer F with optimally tuned 


parameters was applied. It can be shown that these two signals were almost identical in both ” and ¥ 
f \_ £0 /(Z)II. ee 
directions in the sense that | fre(¥) Lc(k)|| is minimized. 


6. CONCLUSION 

In this paper, we have proposed tuning algorithm to improve transparency through impedance 
matching between master and slave manipulators in the multilateral teleoperation system. The passive 
decomposition technique has provided substantial approach to convert the total dynamics of a multi-slave 
system into the bilateral teleoperation architecture so that the existing bilateral control technique can be used. 
Subsequently, from the result of decomposition and the initial states of all slave agents, the teleoperation 
time, 7},.;. that separated the two operations (grasping and teleoperation) has been estimated. In terms of 
obtainining the optimal equalizer, we conclude that FRIT provides powerful algorithms since only a single- 
shot of measurement data was required for tuning. The overall performance and comparison presented 
through numerical example has illustrated the applicability of the proposed tuning method. 
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